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CN ' Abstract 

^P . A mechanism for lifting the cosmological upper bound on the axion 

decay constant, fa, is proposed. It entails the near masslessness of the 
radial mode whose vacuum expectation value is fa- Energy in the 
coherent oscillations of the axion field in the early universe gets fed 
_ , into the motion of the radial mode, from which it is then redshifted 

Qj ' away. It is found that the initial value of fa can be at scales between 

,i^ ■ 2 X lO^'' GeV and 10^^ GeV. This evolves with time to values close to 

the Planck scale. It is suggested that the nearly massless radial mode 
might play the role of quintessence. 
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1 Introduction 

The axion^ idea is an extremely elegant approach to solving the Strong CP 
Problem. However, it is not without difficulties. The main difficulty is that 
axions have neither been observed directly in the laboratory, nor indirectly 
through astrophysical effects. This implies that the axion decay constant, fa, 
must be greater than about 10^° GeV. Moreover, there is the cosmo logical 
bound, based on the so-called "axion energy problem" ,^ that fa is less than 
about 10^^ GeV. This means that the axion decay constant (or, equivalently, 
the scale at which the Peccei-Quinn symmetry is broken) cannot be at any 
of the mass scales at which other physics is known or suspected to be based. 
In particular, it cannot be at the weak scale, the Planck scale, or the grand 
unification scale. 

A number of attempts^ have been made to weaken or remove the cosmo- 
logical bound on fa- In this paper we suggest an idea for doing this. The idea 
is that if the radial mode associated with the axion (that is, the scalar field 
whose expectation value is fa) has a nearly fiat potential then it would have 
increased with time as soon as the coherent axion oscillations commenced. 
As a result of this, energy would have drained out of the axion oscillations 
into the nearly massless radial mode, and from thence been redshifted away 
by the cosmic expansion. We find that the initial value of fa can be as large 
as 10^^ GeV. In the course of cosmic evolution this would have "relaxed" to a 
value we presume to be near the Planck scale. Axions in this scenario would 
therefore be extremely hard if not impossible to detect. However, the radial 
mode may play the role of "quintessence".^ 

The plan of the paper is as follows. In section 2, a rough sketch of the 
axion energy problem in conventional axion models is given. In section 3, 
a simplified discussion of the relaxing axion scenario is given to show how 
energy is drained out of the axion field. In section 4, a more detailed quan- 
titative analysis is given, still assuming, however, that the radial direction is 
exactly flat. In section 5, the radial potential is discussed, in particular how 
fiat it must be and how such flatness might arise. 

2 The usual axion energy problem 

In ordinary axion models, the Peccei-Quinn symmetry is spontaneously bro- 
ken by a complex scalar field, which we shall denote $. The expectation value 



of $ is called fa, while the phase of $ is the axion mode. Thus one may write 
$ = re*^, where (r) = fa, and the axion field a is given by a = faO- Since the 
Peccei-Quinn symmetry is anomalous, QCD instanton effects give the axion 
a mass. The instanton-generated potential has the form Vj = —fiQCOs{N6), 
where /iq is of order the QCD scale, and A^ is an integer, which for present 
purposes can be taken to be 1. 

When the temperature in the early universe was well above the QCD 
scale the instanton-generated potential was not fully turned on, whereas 
after the temperature fell below the QCD scale it had essentially its full 
zero-temperature strength. Let us call the time when the instanton potential 
reached nearly full strength, to. Since the temperature was then of order /io, 
to ~ Mp/fiQ. Assuming a prior epoch of inflation, the axion field at to was 
approximately spatially constant, but barring an extremely unlikely coinci- 
dence, it had no reason to be sitting at the minimum of its potential (which 
we take to be at ^ = 0). Rather, it had some arbitrary value ^(to) = ^q. 
Consequently, the axion field underwent coherent oscillations about its min- 
imum. The energy density in these oscillations at time to would have been 
approximately I/Xq^o- '^^^^ energy scaled as R~^, where R is the scale factor 
of the universe, and therefore the energy density in the coherent axion oscil- 
lations remained (after to) in a constant ratio to the baryon energy density. 
Since at to the baryon energy density was of order 10~^°mp/iQ, the present 
ratio of axion to baryon energy is roughly {^9^(1^) / {lO'^^nipfil) ~ 10^6'q. If 
the axions are not to overdose the universe, 6q ^ 10^^. 

This bound on 6q can be satisfied if fa is sufficiently small. The point 
is that the instanton-generated potential of the axion field did not turn on 
instantaneously. If the potential turned on slowly then the number of quanta 
in the coherent axion oscillations remained constant since it is an adiabatic 
invariant. This implies that as the mass of the axion increased the amplitude 
of the oscillations decreased, and therefore by to, when the potential reached 
full strength, the amplitude 6q could have been very small. A simple estimate 
gives that 9q ~ fa/Mp. Thus, to satisfy the bound on Paxion one requires 
that fa ;S 10^2 Q^y 



3 The relaxing axion mechanism 

Let us suppose that the field $ whose phase is the axion has a fiat potential 
in the radial direction. Then the Lagrangian density can be written 



(1) 

As before, $ = re*^, and the axion field is a = r6. Ignoring spatial derivatives 
of the fields, the equations of motion of r and 6 are 

= r + 3Hr-re^, 

(2) 

= e + {3H + 2^)e + i2lr-^e. 

The last term in the equation for the radial mode r is just the centrifugal 
force, and it is this that drives r to larger values as 6 oscillates. 

The energy densities in the radial and axion modes are given simply by 
Pr = \r'^ and pe = \r'^6'^ + \pq9'^- The energies in a co moving volume in 
these modes are given by E^ = R^pr and Eg = R^pe- Using the equation of 
motion of 9 to eliminate 9, it is easy to show that 



Eg = R 



ir^02 + SHi-^r'9' + ^p',9') 



(3) 



If the oscillator parameters evolve adiabatically, then averaged over many 
oscillations the kinetic and potential energy in the oscillator should be equal. 
That is, (|r^^^) = {\pq9'^) = \{pe)- Therefore, averaged over many oscilla- 
tions, 

Ee = --Eg (4) 

r 

In a similar way, using the equation of motion of r to eliminate f , and aver- 
aging over many oscillations, one finds that 

Er = -Eg - SHEr. (5) 

r 
The interpretation is clear. The increase in r caused by the centrifugal force 
drains energy out of the axion oscillations and into the radial mode, while 
at the same time energy is redshifted away from the radial mode. To put it 
another way, as r increases the effective mass of the axions, Pq/t, decreases, 
while the number of axions remains constant. Thus Eg ~ 1/r, as implied 
also by Eq. 4. 



Suppose that r ~ f^, and Eg ~ t^^. The exponent q can be determined by 
writing the equation of motion of r in terms of Eg as follows: = f + 3Hr — 
R^^{EQ)r^^. The first two terms go as t^''~'^\ while the last goes as t^('^i+^/'^) ^ 
assuming that R ~ t^/^. Therefore q = 1/6. Writing Eg = EgQit/to)"^/^, 
Eq. 5 is solved by E^. = ^Egoit/to)'^^^ + ct~^/^. Thus, for large t, the radial 
energy is one-eighth of the energy in the coherent axion oscillations, and the 
energy in a comoving volume in either mode falls off as t~^^^. 



4 A more detailed analysis 

So far we have not taken into account the temperature dependence of the 
instanton-generated potential of the axion field. According to the dilute- 
instanton-gas calculation of Gross, Pisarski, and Yaffe,^ the instanton po- 
tential for 6 at high temperature goes as /i^ ~ T^ exp(— 87r^/5'^(T)) ~ T^ 
exp(i(llA^ - 2Nf)\nT) ~ T^-'^+^^f^^l We shall assume henceforth that 
H^ = ^liTo/Tf' = /i^(t/to)^ for T > /io (i.e. t < to), and that /i^ ^ /i^ for 
t ^ to- The dilute-instanton-gas calculation suggests that k = 5/2, but we 
shall keep /c as a parameter. 

If we change variables to r = t/to, and denote d/dr by a dot, then we 
can write the equations of motion 

= r + ^r - r(P-, 

At ' 

= ^- + (I: + 2i)^ + (/i^tg/r V^. (^) 

The effective mass of the field ^, then, is mg{T) = (/iQto/'r(r))r'^/^. Therefore 
it is reasonable to make the ansatz that 6{t) has the form 

0(r) = A(r)e^^M, (7) 

where A and B are real and 

Bir) = mg{r) = (^^to/r(r))r'=/l (8) 

Substituting into the equation of motion for 6 and taking the real and imag- 
inary parts of the equation, one gets (using the fact that r/r = k/2T — B/B) 



A 



k)/T - 2B/B 



= -B+ U + k)/T + 2A/A 



A, 
B. 



(9) 



The first of tliese equations is exactly solved by (A/B'^)t^^^'^^''^ = constant, 
while the second is solved exactly by (B/A'^)t~'^^^'^~^''^ = constant. Eliminat- 
ing B from these equations and solving gives 



where c', C, and Oq are integration constants. After many oscillations the 
integration constant C can be neglected, and one has, using Eqs. 7 and 10, 
a solution for 9{t) of the following form: 

e(T) = ^or-(5/6+'=/3) exp(^6or ('/'+'=/') + i^). (11) 

(^0 = q/R I fr/s O Elqs. 8 and 10 directly give the solution for the radial mode: 



5/6+k/3^' 
(5/6+fc/3)bo 



r^r) = (5/6^°fc/3)j) r*^"^/^^^'/^^ (again, neglecting the integration constant C). 
Note that for the case of A; = 0, which corresponds to a fixed value of /x^, 
the radial variable increases as r^^^, in agreement with the result obtained 
in section 3. 

It must be checked that this solution for r(T) and the solution for 6(t) 
given in Eq. 11 satisfy the equation of motion of r. Assuming that 6 is 
dominated by the rapid oscillations of 6 rather than by the slow variation of 
its amplitude, one has that 9 = ibo{5/6 + k/3)eoT-^ exp{iboT^^/^+''/^'> + ih'o). 
Averaged over many oscillations, therefore, {6"^) = |6q(5/6 + k/SyO^r'^^. 
Substituting this and the expression for r{T) into the equation of motion for 

r (see Eq. 6), one sees that that equation is satisfied if b^ = ^ ^ ^^ 9^^ . 

Thus we have that 

T{r) = , ^^°^° e,r^'/'+'/'\ (12) 



,/0^^kj{ATk) 



The solutions given in Eqs. 11 and 12 apply to the period r < 1, in which 
the instanton potential was still turning on. The same expressions with k 
set to zero apply to the period r > 1, after the instanton potential turned 
on. The true solution will smoothly interpolate between these in the period 
r ~ 1. (Note that the k ^ and A; = solutions for 6 actually agree at 
r = 1, while the k ^ and A; = solutions for r differ at r = 1 by a factor 



of y'(l + A;)(4 + A;)/2^2.4.) 

One is now able to estimate the energy in the axion oscillations. The 
crucial parameter is the value of r at the time when the axion oscillations 
started. We will call that time tj. The oscillations started when the effective 
mass of the 6 field, rriQ = B, became equal to the expansion rate of the 
universe, H = (2t)-^ That is, when {2ti)-^ ^ {l4'tohiU)){U/tof'^- Clearly, 
the smaller r(tj) was, the earlier the axion oscillations began. Turning this 
around, r(t,) ^ 2/i2to(t,/to)(i+^/2) ~ Mp{U/toY^+^'^\ 

At ti it is to be expected that 9 was of order unity, since it had not had 
time to be affected by the instanton potential. But according to Eq. 11, 
0{ti) ~ eQ{ti/tQ)-^^/^+^l^\ Thus, ^0 is of order {ti/tQY^'^+'''^\ or, in terms of 



9o ~ (r(t,)/Mp)(^+^). (13) 

As was seen in section 2, the factor 6q tells how much the energy of the 
coherent axion oscillations was suppressed by the time the axion potential 
fully turned on at to- We will call this suppression factor Sbefore- 



Sbefare ~ 0^ ~ (r (t,)/Mp) 3 ^ Vf). (M) 

If this were the only suppression of the axion energy, solving the axion energy 
problem would require that Oq ;$ 10^^^^. With k = 5/2 this gives r(tj) < 
2 X 10^^ GeV. That is, the initial value of "/„" can have been quite near the 
grand unification scale. By to this would have increased, according to Eq. 
12, to a value r(to) ~ /io^o^o ~ 9qMp ~ 3 x 10^^ GeV. One possibility, which 
we will call Case I, is that the radial mode stopped evolving at that point, 
because its potential has a minimum there. Another possibility, which we will 
call Case II, is that r continued to increase to some final value near the Planck 
scale. That would mean that the axion energy would have been further 
suppressed by the evolution of r in the period t > to. Since r(to) ~ OqMp, 
and it is assumed that r{tf) ~ Mp, there is an increase of r by a factor of 
9q^ in this period. It is easily shown that the energy of the axion field in a 
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comoving volume varies inversely with r in this period (as was seen already 
in section 3), so that the further suppression of the axion energy, which we 
shall call S after i is given by 

S after ~ ^0) (15) 



or 



Jtotal — JbeforeJafter ~ "o- \^^) 

In Case II, therefore, it is only necessary that 9q < 10^^/^, meaning that 
r(tj) < 10^^ GeV. Case I and Case II are, in a sense, the extreme cases. One 
can consider intermediate cases as well. But one sees that in general the 
relaxing axion scenario would have fa starting out in the range 10^^ to 10^^ 
GeV, near the grand unification scale, and evolving to higher values. 

5 The flatness of the radial potential 

The mechanism described in the preceding sections depends crucially on the 
assumption that the potential in the radial direction is nearly flat. For the 
mechanism to work, the centrifugal term in the equation of motion for r had 
to have dominated over the force coming from the potential energy of r. That 
is, r9^ > \V'{r)\. One can write this as paxion > \rV'{r)\, where paxion is the 
energy in the coherent axion oscillations. 

At some point this condition was no longer satisfied and the radial field's 
evolution was controlled by V{r). Unless the radial field was at that point 
overdamped (not so in the cases of interest) it would have started to oscillate 
about the minimum of V{r). For reasonable potentials (where one assumes 
that the cosmological constant problem has somehow been solved) one would 
expect that V{r) ~ |rV^'(r)|, and therefore the energy in the coherent oscil- 
lations of the radial mode were also of that order. Consequently, when the 
coherent radial oscillations began, the energy in them was typically of the 
same order as the energy in the coherent axion oscillations. Thus the coherent 
radial oscillations do not in themselves pose a cosmological problem. 

However, for the mechanism to work at all it is necessary that the cen- 
trifugal term in the equation of motion of r dominated over the potential 
term for a sufficiently long time to solve the axion energy problem. This 
puts a constraint on the flatness of V{r). 



In Case I, it is assumed that the centrifugal term drove r until to, when 
the temperature was of order /iq. At that time paxion had to have been less 
than about lO^p^ ~ IQ^^ji^mp ~ lO^^'^GeV. Thus, at to it must also have 
been that \rV'{r)\ < IQ-^^GeV^ 

In Case II, the centrifugal term is assumed to have dominated until r got 
to be of order Mp. Since r(to) ~ OqMp, and r grew as t^/^ for t > to, this 
happened at a time tj ~ ^(7^^o- Assuming a radiation dominated universe, 
T{tf) ~ ^QT(to) ~ StotaiPo ~ 10~^/io- Therefore, at t/ it must have been that 
\rV'{r)\ < lO-^T{tffmp ~ IQ-^^GcV^ 

How can the potential be that flat? Certainly it is trivial to arrange that 
the radial direction be flat in the supersymmetric limit. The real problem 
is to insulate the radial mode from supersymmetry breaking. This is easiest 
to do if supersymmetry is broken at low energies, as it is in models with 
gauge-mediated supersymmetry breaking.^ In such a model, the ordinary 
quarks must be split from their supersymmetry partners by an amount that 
is of order 10^ to 10^ GeV, which scale we will call niQ. Since the axion 
sector must couple directly or indirectly to the quark sector for the Peccei- 
Quinn symmetry to have a QCD anomaly, supersymmetry breaking will be 
fed into the radial mode of the axion sector through quark loops. Typically, 
then, the radial mode of the axion sector will acquire a potential of the form 
e?72Q ln(r/mo)- e is model-dependent and is smaller the more indirectly and 
the more weakly the axion sector couples to the quark sector. In Case I, one 
has that em^ < 10^^°GeV^ implying (if mo ~ 1 TeV) that e < 10~22_ j^^ Q^gg 
II, one has the more severe constraint that e ;S 10^"^^. 

To see how such small values of e might be achieved, consider first a 
conventional axion model where the radial mode has a tree-level potential. 
The relevant terms in the superpotential would have the following general 
structure: 

Wax^on = QSQQ + Ws, (17) 

where Q and Q are lefthanded quark and anti-quark superfields, and 5 is a 
superfield containing the axion. Ws is some set of terms, generally involving 
other fields, which has the effect of fixing | (5") | to have some value M. The 
phase of S is, however, assumed not to be fixed except by QCD instanton 
effects. One can write S = (M + S')e*^, where here we mean by S the bosonic 
component, and the axion field is M6. The radial mode S would typically 
have some mass of order M. 



Consider, now, a somewhat different model, with the corresponding terms 
in the superpotential being the following 

Wa.,on = gSQQ + Ws + g'iSA - M'B)Y, (18) 

where S, A, B, and Y are all gauge singlets, and as before Ws has the effect 
of making \{S)\ = M but leaving the phase of S undetermined. Suppose 
that both M and M' are very large compared to the scale of supersymmetry 
breaking. It is apparent that, since S has a Peccei-Quinn charge, so must 
either A or B. Therefore, if (A) and (B) are larger than M, they rather than 
{S) control the value of fa- 

Assume that Y has no other couplings in the superpotential. Then one 
of the terms in the scalar potential is \g'{SA — M'B)\ . This term fixes only 
the ratio of A and B and leaves them otherwise undetermined. There is, 
therefore, a direction that is flat in the supersymmetric limit along which A 
and B can become arbitrarily large. This flat direction would play the role 
of r in this model. 

Writing S = {M + S)e'^, A = re*", and B = (^r + 5)e^("+''+''), the 
aforementioned term can easily be found to give a mass to the fields S, B, 
and 6 = r6. In particular, one finds that \g'{SA — M'B)f = g'^M'^P + 
g '^{{r)S — M'B)^+ terms of cubic and higher order. In addition, there is the 
mass term for S coming from Ws- The phase a is an exact goldstone mode, 
while 9, since it corresponds to an anomalous U{1) by virtue of its coupling to 
the quarks, is the axion mode. The scalar field r has (before supersymmetry 
breaking) a flat potential. It is easily seen that the Lagrangian terms for 
r and 6 have (after suitable rescaling of fields) essentially the same form 
as those shown in Eq. 1. Thus this model is an implementation of the 
relaxing axion mechanism. The question is how large a potential r gets in 
this model. If supersymmetry breaking is mediated from some hidden sector 
by Standard Model gauge interactions, then only Q and Q will directly feel 
it. The splittings in the supermultiplet S will only arise through a quark 
loop, and the splittings in the multiplets A and B will only arise through 
diagrams involving both an S loop and a quark loop. One would expect 
that e ~ ((7^/167r^)((y' ^/IGtt^). This can easily be as small as required. For 
example, if M and M' are near the unification scale, the coupling g could 
even be as small as 10~^^. 

It seems to be possible to shield the radial mode from supersymmtry 
breaking even more thoroughly. As an extreme case one could imagine a 
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whole series of gauge-singlet sectors separating the quark sector from the 

axion sector: Squark < ^ 'S'l < > ... < > Sn < ^ Saxion, where the S"s denote 

various sectors, and the arrows represent a coupling in the superpotential 
between two sectors. In this case, supersymmetry breaking in the radial 
mode would be an n-loop effect. On the other hand, because of the couplings 
between sectors, there are fields in each sector that have non-trivial Peccei- 
Quinn charges. Thus the axion field, though removed by several steps from 
the quarks, will nonetheless get a potential from QCD instanton effects that 
goes as — /ig cos{N6), where A^ is an integer. 

From these examples, it seems that there is no reason in principle why 
the radial mode could not be sufficiently fiat to allow the relaxing axion 
mechanism to work. It is an interesting question whether this flat radial 
direction can be identified with other fiat directions that have been discussed 
in recent years. Could it be, to mention two obvious examples, the dilaton 
of superstring theory or the "quintessence" field? In any event, given the 
increasing role that very fiat potentials have played in particle physics and 
cosmology (infiatons, quintessence, dilatons, moduli, and the axion itself), 
it is suggestive that a fiat radial direction can allow the breaking of the 
Peccei-Quinn symmetry to take place at the unification scale or even higher. 
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